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Electrical modulation of the edge channel transport in topological insulators 
coupled to ferromagnetic leads 

Yuan U, 1 ' 2 ^ M. B. A. Jalil, 1 ' 3 Seng Ghee Tan, 3 and GuangHui Zhou 4 

^Department of Electronic and Computer Engineering, Information Storage Materials Laboratory, 
National University of Singapore, 1 Engineering Drive 3, Singapore 111516, 



^Department of Physics, Hangzhou Dianzi University, Hangzhou 310018, P. R. China 

^ Data Storage Institute, DSI Building, 5 Engineering Drive 1 (Off Kent Ridge Crescent), 

National University of Singapore, Singapore 111608, Singapore 

^Department of Physics and Key Laboratory for Low- Dimensional Quantum Structures 
and Manipulation (Ministry of Education), Hunan Normal University, Changsha 410081, 
China 

(Dated: 16 August 2012) 

The counterpropagating edge states of a two-dimensional topological insulator (TI) carry electrons of opposite 
spins. We investigate the transport properties of edge states in a two-dimensional TI which is contacted 
to ferromagnetic leads. The application of a side-gate voltage induces a constriction or quantum point 
contact (QPC) which couples the two edge channels. The transport properties of the system is calculated 
via the Keldysh nonequilibrium Green's function method. We found that inter-edge spin-flip coupling can 
significantly enhance (suppress) the charge current when the magnetization of the leads are anti-parallel 
(parallel) to one another. On the other hand, spin-conserving inter-edge coupling generally reduces the current 
by backscattering regardless of the magnetization configuration. The charge current and the conductance as 
a function of the bias voltage, also exhibit similar trends with respect to spin-flip coupling strength, for both 
parallel and anti-parallel configurations. Hence, gate voltage modulation of edge states via a QPC can provide 
a means of modulating the spin or charge current flow in Tl-based spintronics devices. 

PACS numbers: 73.43.-f, 72.25.Dc, 85.75.-d 
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I. INTRODUCTION 

Topological insulators (TIs) are electronic materials 
that possess an insulating bulk gap and topologically- 
protected conducting states on their edges or surfaces^. 
Bi2Se3 and related materials with gapless surface 
states*—, e.g., Bi2Te3 and Sb2Te3, have recently at- 
tracted extensive attention in the condensed matter and 
device physics communities^—. Thin TI films were pro- 
posed for application in the magnetic memory-^ due to 
their extremely high surface-to-volume ratio^i. Accord- 
ingly, a mechanical exfoliation method was developed 
for preparation of thin TI films with significant surface 
transport^—. 

On the other hand, the two-dimensional version of 
topological insulators is also known as a quantum spin 
hall insulator, consisting of edge states that behave 
as perfect one-dimensional counterpropagating (helical) 
channels. These channels are well-localized along the 
edges, carry electrons with opposite spins, and are pro- 
tected from nonmagnetic impurity scattering by time- 
reversal symmetry. The peculiar properties of the he- 
lical edge states have been demonstrated theoretically 
and experimentally in HgTe/CdTe quantum wellsi^— , 
and have been predicted to exist in various other 
materials^—. Due to the robust spin property of these 
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edge channels, two-dimensional topological insulators are 
promising candidates for spintronic devices. 

In any future spintronic device, it is desirable to 
achieve electrical (rather than magnetic field) manipula- 
tion of spin transport. However, the robustness of the he- 
lical edge-state transport to nonmagnetic impurity scat- 
tering and local perturbation, means that it is not easy 
to electrically modulate the edge-state transport in TI- 
based devices. Nonetheless, investigations into electri- 
cal control of edge channel transport in TIs have gath- 
ered pace recently. For instance, due to the finite size 
effect^, the inter-edge tunneling arising from the over- 
lap between states from opposite edges can open up an 
energy gap, and modify the conductance of the system. 
The amplitude of the inter-edge tunneling is determined 
by the finite decay length of the helical edge states into 
the bulk>2£r— . This phenomenon may provide a means 
to adjust the transport property of edge states, via finite 
size effect of the quantum point contact^, interferom- 
etry of the edge states^, and nanoscale engineering of 
the edge geometry, e.g., by patterning a constriction in a 
HgTe heterostructure^ 3 .. 

In this article we propose the design of a Tl-based 
spintronics device which incorporates a two-dimensional 
topological insulator (2D-TI), e.g. a HgTe/CdTe quan- 
tum well, contacted to two ferromagnetic leads. We 
shall focus on the electrical manipulation of the trans- 
port properties by inducing inter-edge tunneling through 
an electrically induced quantum point contact. This can 
be achieved by applying external voltage via a pair of 
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FIG. 1. Schematic description of the proposed two-terminal 
setup, where edge states flow at the top and bottom bound- 
aries of the topological insulator. There exists a QPC in the 
central region which can be controlled by the split-gate volt- 
age V g . Black (red) lines refer to spin- up (spin-down) edge 
channels, which is assumed to be polarized along the z di- 
rection. Application of the split gate voltage V a reduces the 
spatial separation of the two edge states, thus increasing the 
inter-edge tunnel coupling at the QPC region, which is de- 
picted by the blue (green) dashed lines. 



split gate along the transverse direction of the device. 
We apply the nonequilibrium Keldysh Green's function 
method to calculate the charge current flowing through 
the Tl-based device. It is found that conductance of such 
device can be significantly enhanced by the presence of 
inter-edge spin-flip tunneling when the leads are highly 
spin asymmetric and their magnetizations are in the an- 
tiparallel configuration. Thus, the proposed setup can 
potentially realize an electrically controlled TI device in 
which the conductance varies according to the inter-edge 
tunneling coupling and the magnetization configurations 
of the leads. 

The organization of the rest of the paper is as follows. 
In Sec. HU the Hamiltonian of the system is introduced, 
and the current formula is derived based on the nonequi- 
librium Green's function method. In Sec. IIII1 we present 
the results of our numerical calculation of the charge cur- 
rent and the differential conductance in the presence of 
inter-edge tunneling couplings modulated by the split- 
gate voltage. Finally, a brief summary is given in Sec. 
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II. MODEL AND FORMULA 



We consider a Tl-based device consisting of a channel 
made of a 2D-TI bar sandwiched between two ferromag- 
netic leads [see Fig. [T]. We assume the presence of a 
gate-bias induced QPC, across which inter-edge coupling 
occurs. In this paper, we use the low energy effective 
theory of the edge states to investigate the transport 
property of the QSH syste m 32 i 34 . This system can be 
described by the following Hamiltonian: 



where He is the effective Hamiltonian describing the 
QSH edge states 

He = J2(vpV^k)cl k<T Cfj ka + Yy v cic\ ka c bka 

/3fcer key 

+Vv V c2c\ k(J Cbka + C.C.]. (2) 

In the above, Cg fccr (cpku ) is the creation (annihilation) 
operator of spin a =t (4-) for electron constrained along 
the /3 edges, where the edge index (3 = t,b refers to the 
top- and bottom-edge channels. The first term in He 
describes the left- and right-moving QSH edge states, 
while the second and third terms represent, respectively, 
the spin-conserving and spin-flip couplings between the 
two edge channels. Accordingly, v is the edge state ve- 
locity, V c i and V C 2 are, respectively, the spin-conserving 
and spin-flip tunnel coupling strengths between the edge 
states. These can be modulated by means of the split- 
gate voltage V g which modifies the width of the QPC 32 . 
The edge and spin states are characterized by the fol- 
lowing indexes: r\ t j b — ±1 and f\\n = ±1. In general, 
in order to preserve time-reverse symmetry, only the two 
types of coupling represented by V c i and V C 2 are allowed 2 . 
In the device operation, an electron is injected from the 
source electrode and flows through along one of the edges 
of the topological insulator. An important point to note 
is that under nonequilibrium condition, the total current 
is contributed by electrons with different wavevector k 
corresponding to the conduction "window" between the 
electrochemical potentials of the source and drain elec- 
trodes. However, for simplicity, one can assume the cou- 
pling constants to be independent of k and energy over 
the conduction window (this is especially valid for small 
source- drain bias). 

The Hamiltonian Hl and Hr, which describe the left 
and right ferromagnetic leads, and Ht which describes 
the coupling between the FM leads and the central TI 
region are given by: 
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{(-ka + V&M) is the energy of conduction electrons in the 
a lead, and is characterized by the amplitude k of the 
wave vector, and the amplitude M and orientation r\ a of 
the magnetization of the FM lead. a kaer and akau denote 
the creation and annihilation operators of electrons with 
spin a in lead a, while Vk a ,fia denotes the tunnel coupling 
between the edge state of the TI and the lead. 
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The charge current can be calculated using the stan- 
dard Keldysh nonequilibrium Green's function method 3 -^, 
namely, 



Jot 



de 
2^ 



Tr{r Q [G<( £ ) + U(e)(G r (e) G a ( £ ))]}, r 



(5) 



where / Q (e) = {exp[(e — ^ a )/kBT] + is the Fermi- 
distribution function, and the trace is over both the spin 
and edge-state indexes. The linewidth r Q (e) is repre- 
sented by a (4 x 4) matrix in spin space and edge-state 
space, which can be written as 



r a (e) 



with 




r' 3 o 

1 at 



u a , 
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u a = 



cos(6> Q /2) sin(0 a /2) 
-sin(0 a /2) cos(0 Q /2) 



Here, 8 a denotes the magnetization direction in lead a, 
so that we have 8 L = and 8 R = 8, i.e., fixed (variable) 
magnetization orientation for lead a = L (R). The com- 
ponents of the linewidth matrix are related to the tun- 
nel couplings as follows: T^ a (e) = 2irJ2 k \Vk a ,/3a\ 2 S{e - 
£kaa)- In the following numerical calculation, we assume 



am& = Kt/AD - |e|) and 1*.^ = ^T (l±p a ), 
where D = 500 is the bandwidth, T = 1 is the unit of 
energy, p a denotes the spin asymmetry factor for tunnel- 
ing across the barrier, and £ a is the asymmetry in the 
Tl-lead coupling between the left and right barriers. As- 
suming the junctions are symmetrical and the leads are 
made of the same FM material, we have £l = £,R. = 1 
and Pl — Pr — P- The retarded Green's function G r (e) 
in the current formula of Eq. ([5]) is also a (4 x 4) matrix, 
which can be solved via the standard equation of motion 
technique 3 ^, which yield the following Dyson's equation 

G r (e) = G r '°(e) + G r >°(e)(£ r ( L ) + 17^ + V c )G r (e). 

(6) 

G r '°(e) describes the edge states of the topological insu- 
lator in the absence of TFlead and inter-edge tunneling 
couplings, i.e., 



G r '°(e) = 



1 



(e - vk Q i]y + ii]) 








(e + vk rjv + it]) 








(e + vko7}v + irf) 








1 



(7) 



(e - vk riv + irf) I 



where kg denotes the wavevector associated with the 
Fermi energy, V = (hl — Hr.) is the bias voltage and 
rjv = ±1 is adopted to describe the helical property of 
the edge channels. £ r / a ( Q ' = (^i/2)T a (e) are the self- 
energies due to the lead a, while the matrix V c denotes 
the self-energy due to inter-edge tunneling, and is explic- 
itly given by 



/ 



v r „ = 







V cl V c2 \ 



-v c2 v cl 

V c i -V c2 



\v c2 v cl 
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Finally, the lesser Green's function G*- in Eq. ([5]) can 
be calculated using the Keldysh equation: G < (e) = 
G r (e)(£ <(i) + £ <(i?) )G a (e), with the scattering func- 
tion being given by £<( Q ) = if a T a . After solving for 
the charge current J a , the differential conductance Gd 
can be calculated straightly, Gd = dJ/dV with Jl = 



—Jr = J. Accordingly, the tunneling magnetoresistance 
(TMR) can be derived by the ratio (G d (O)-G d (0))/G d (6) 
with G d (8) denoting the conductance associated with the 
angle 8 of the right lead's magnetization. 



III. NUMERICAL CALCULATION 

In this section, we numerically investigate the trans- 
port property of the Tl-based device. We adopt the 
parameters of HgTe quantum wells given in Ref. @, 
i.e., with thickness d = 7nm, A = 364.5meV nm, B = 
— 686meV nm 2 , M — — lOmeV, and D = 512meV nm 2 , 
while the edge state velocity is v ~ 5.5 x 10 5 m/s. The 
width of the device is set at W = lOOOnm so that 
the energy gap of edge states is very tiny and can be 
neglected 2 ^. The pinching of the two edge states re- 
sults in a local modification of the spin-orbit coupling, 
which may induce spin-conserving and spin-flip tunnel- 
ing at this narrow regio n 4 i 32 i 33 . We adopt the two phe- 
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FIG. 2. (a) and (b) Charge current J as a function of the 
angle 9 under different inter-edge tunnel coupling conditions. 
In (a), we analyze the role of the spin- flip inter-edge cou- 
pling, characterized by V C 2, while in (b) we consider the effect 
of spin-conserving inter-edge coupling, characterized by V c i. 
Other parameters are To = 1, T = 0.2, p — 0.8, [io = 0.3 and 
V = 0.5. 



nomenological parameters V c \ and V C 2 to describe the two 
tunneling processes and their relative amplitudes. The 
Dirac point at equilibrium is set at energy E = 0, while 
the Fermi energy is at /i = vk = 0.3r (r = lOmeV). 
A bias voltage V is applied to the device, such that the 
chemical potentials of the leads are fiL — Mo + V/2 and 
V-R = Mo - V/2. 

We first analyze the band dispersion of the edge states 
in the presence of the inter-edge tunneling couplings. 
From the Hamiltonian He of a topological insulator, the 
energies of top-edge states are given by E a — ±vk in 
the absence of inter-edge tunneling couplings. Once the 
split-gate voltage V g is applied and inter-edge tunnel cou- 
pling becomes finite, the energies of the edge states are 
modified to 

E' a = ±yJ(vk + V c2 )* + V*. (9) 

The energies of the bottom-edge states can be calculated 
using the same method. The spin-preserving part of the 
inter-edge tunnel couplings causes a band gap opening 
of A = 2V c x. Thus electrons cannot flow through the 
device if the Fermi energy is adjusted to lie within the 
energy range of |mo < V c \. In this paper, however, we 
focus on the transport property of the Tl-based device in 



FIG. 3. Bias dependence of the charge current for the cases (a) 
8 — 0, (b) 6 = 7T. The upper insets in (a) and (b) show the bias 
dependence of the differential conductance, while the lower 
inset in (b) is the corresponding plot of the TMR. V c i = 
and the other parameters are the same as those of Fig. [2] 



the regime of |mo| > V c i, where there exists current flow. 
We initially study the effect of the inter-edge tunnel cou- 
pling on the charge current. Fig. Ufa) shows the effect 
of the spin-flip tunneling, characterized by the coefficient 
V C 2 on the charge current as a function of magnetiza- 
tion orientation 9. The charge current decreases with 
increasing tunneling strength V C 2 when 9 < 5tt/12, es- 
pecially at 9 = [the parallel configuration]. However, 
there occurs an opposite trend, i.e., the charge current 
can increase with increasing V C 2 when 9 > 5tt/12. This is 
especially so when 9 = ir [the anti-parallel configuration] . 
The observed trend can be explained by considering the 
spin-state of electrons in the edge states. In the parallel 
configuration, spin polarized electrons can flow through 
the central TI region along one of the edge channels in 
the absence of the inter-edge coupling. For instance, in 
the schematic diagram of Fig. 1, we find that when the 
magnetization of the leads are pointing in the z-direction, 
then most of the forward electron flux will flow in the top 
edge state. In the presence of inter-edge spin-flip cou- 
pling, however, some of the forward moving electrons in 
the upper edge will be scattered to the spin-down state 
in the lower edge. They will encounter a higher probabil- 
ity of being blocked at the right electrode since they are 
a minority carrier there. On the contrary, for the anti- 
parallel case, the current is suppressed in the absence of 
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the inter-edge coupling. Referring to Figure 1 again, let 
us assume that the magnetization of the left (right) leads 
are magnetized in the +z (—z) direction. The forward 
spin-up flux in the top edge will be blocked at the right 
electrode as spin-up electrons are minority carriers there. 
The presence of inter-edge spin- flip coupling (i.e., finite 
V c i) enables the spin- up electron in the top edge to be 
scattered into the forward moving spin-down state in the 
lower edge. The spin-down electron will have a higher 
transmission into the right electrode as it will be a ma- 
jority carrier in that electrode. Thus, by utilizing the 
presence of spin-flip tunneling at the QPC, one can mod- 
ulate the transport properties of this Tl-based device for 
different magnetization configurations. 

At the same time, the spin-conserving inter-edge tun- 
neling coupling at the QPC induces the back-scattering 
process. Since this can coexist with the spin-flip tun- 
neling term, it is deservable to investigate the combined 
effect of these two terms on the charge current, as shown 
in Fig. EJb) . When there exists only the spin-conserving 
tunneling, for example V c i = 0.3 and V C 2 = [red dashed 
line], the charge current is decreased by the tunneling 
term via the back-scattering process. In other words, 
some of the forward moving electrons in the upper edge 
will be scattered to the backward moving spin-up states 
in the lower edge, thus decreasing the overall charge cur- 
rent. The reduction in current due to the back-scattering 
is greater for the parallel configuration than that for the 
antiparallel case. This is because the current reduction 
in the parallel case is mainly due to the back-scattering 
of majority spin carriers, and this involves more electrons 
than that of minority electrons in the anti-parallel case. 

In the presence of both types of inter-edge tunnel 
coupling, i.e. both the spin-flip and spin-conserving 
types, the net result will be suppression (enhancement) 
of charge current when the device is in the parallel (an- 
tiparallel) configuration. In the parallel configuration, 
both the spin-flip and spin-conserving tunneling reduce 
the current, the former by scattering to the minority spin 
state, while the latter by backscattering. However, in the 
antiparallel configuration, the two inter-edge couplings 
have opposite effects — the spin-flip (spin-conserving) 
process increases (decreases) the current. When the two 
are of equal strength, e.g., V c \ = V C 2 = 0.3, the spin-flip 
process dominates, and so there is a net increase in cur- 
rent in the antiparallel configuration [see Fig. HJb)]. The 
current suppression due to the spin-conserving tunnel- 
ing also results in the cross-over between net enhance- 
ment and net suppression of charge current occurring 
at a large-angle configuration beyond = 57r/I2 [see 
Fig.[2fb)]. Thus, the transport properties of our Tl-based 
device can be effectively modulated by adjusting the rel- 
ative strengths of the two types of inter-edge tunneling 
couplings across the QPC. In order to experimentally ver- 
ify the above results, one must first prepare the HgTc- 
based device with transverse width of at least W = 1000 
nm, as schematically shown in Fig. [TJ For the two ferro- 
magnetic leads, the magnetization of the left lead is fixed 



and pointed in the out-of-plane z direction, while the 
right lead is made of a softer magnetic material, so that 
its magnetization direction can be changed from 9 = 
to 8 = 7r with respect to the left lead, by application 
of a magnetic field. One can modulate the width of the 
QPC by using the gate voltage V g , which would in turn 
induce different tunneling strengths V c \ and V C 2, due to 
the finite overlap between edge states on the upper and 
the lower sides^. By analyzing the voltage dependence 
of the charge current, one may determine the variation 
of the charge current as a function of V c \ and V C 2 ■ 

Next, we analyze the charge current J and the con- 
ductance Gd of the Tl-based device under varying bias 
voltage. In Fig. (3{a) and (b), the bias dependence of J 
is plotted for different spin- flip coupling strengths V C 2, 
for the parallel (8 = 0), and antiparallel (8 = tt) configu- 
rations. As noted earlier, for the parallel (antiparallel) 
configuration, J decreases (increases) as V C 2 increases. 
Similarly, the bias dependence plot of Gd also shows 
a decrease (increase) with increasing V c i when 8 = 
(8 = tt), as shown in the upper insets of Fig. 3(a) [(b)]. 
This may be explained by the fact that a finite V C 2 would 
lead to coupling of two edge states with opposite spins, 
resulting in spin mixing and thus reducing the degree of 
spin asymmetry in the device. Consequently, the TMR 
ratio decreases with increasing V C 2, as shown by the 
lower inset of Fig. [2Jb). For both configurations, the 
conductance also decreases with increasing bias voltage 
V, with a larger decrease occurring for the parallel case. 
This translates to a slight decrease in TMR ratio with 
bias voltage V for small coupling strength V C 2- Hence, 
the conductance and the TMR of the Tl-based device 
can be modulated by both the source-drain voltage V, 
and the split-gate voltage V g , which alters the coupling 
strength V C 2 ■ 



IV. CONCLUSION 

In summary, we have investigated the transport of a 
two-dimensional topological insulator (TI) sandwiched 
between two ferromagnetic (FM) leads, and incorpo- 
rating a quantum point contact (QPC) which couples 
the edge states along the top and bottom boundaries. 
The conductance across the device is calculated via the 
nonequilibrium Green's function (NEGF) method. Wc 
found that spin-flip tunnel coupling between the edge 
states mediated by the QPC can enhance (suppress) the 
current when the magnetization of the two FM leads 
is in the antiparallel (parallel) configuration. On the 
other hand, the spin-conserving tunnel coupling leads to 
backscattering which reduces the overall current regard- 
less of the magnetization alignment of the FM leads. The 
conductance and TMR of the Tl-based device can be con- 
trolled by both the source-drain voltage V, as well as the 
split-gate voltage V g , the latter modifying the inter-edge 
coupling strengths across the QPC. Our proposed device 
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suggests a new class of Tl-bascd devices which incorpo- 
rate a QPC to electrically control the inter-edge coupling, 
so as to harness the topological and spin-polarized nature 
of the TI edge states for practical application. 
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